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Abstract 

The production of cosmic-ray electrons of characteristically low energies 
is investigated. Secondary sources, other than that of meson decay, are con- 
sidered, and constraints are placed on both secondary and primary sources. 

(1) ( Calculations are made of the intensity of low-energy knock-on and beta- 

X 

decay electrons which are secondary to cosmic-ray interactions. ' In particular, 
knock-on production is calculated in the 100 KeV to 50 BeV kinetic -energy 
interval. Interstellar losses due to ionization, leakage from the galaxy, synchro- 
tron, bremsstrahlung, and inverse compton effects are considered, as well as 
those due to plasma excitation, the red-shift, bremsstrahlung, synchrotron, 
and inverse compton effects in the intergalactic medium.j It is shown that the 
intensity of low-energy relativistic electrons from these sources is not negligible 
compared with that from the rr-ji-e process, but does not account for the ob- 
served interplanetary electron intensity. 

(2) ^Energy inputs to the injected secondary electrons by a possible solar 
electric field of low magnitude, and by a possible galactic Fermi acceleration 

*NAS-NASA Post-Doctoral Resident Research Associate. 

iii 



are investigated. It is shown that one such input is necessary if the observed 
low-energy interplanetary electron intensity is to be attributed to secondary 
production./ A heliocentric field which does allow for a fit to the low-energy data 
cannot, however, account for the high-energy BeV electrons found to be in excess 
of those from TT-/j.-e production. The Fermi acceleration shown to be necessary 
to provide a fit is greater than that usually postulated for cosmic-ray protons, 
and also requires that the ratio of escape losses to acceleration, k/a , be much 
smaller than is usually assumed for protons. This distinction is acceptable only 
if one postulates a significant difference between interstellar proton and electron 
propagation. 

(3) The observation that the velocity spectrum of electrons in the energy 
per-unit-mass region of 7 to 25 closely approximates that of the cosmic-ray 
protons, and the necessity of constraints on the secondary electron hypothesis 
outlined above, suggest that most of the low-energy electrons are of primary 
origin. The similarity between this conclusion and the conclusion (based on the 
measurement of the charge ratio of electrons) that the higher-energy electrons 
are mostly primary is discussed. 
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ON THE PRODUCTION OF 
LOW-ENERGY COSMIC-RAY ELECTRONS 


I. INTRODUCTION AND DISCUSSION OF OBSERVATIONS 


The study of cosmic radiation has been, for the most part, the measurement 
of the intensities and energy spectra of the protons and other nuclei which pos- 
sess the bulk of the cosmic-ray energy content. Recently, the electromagnetic 
component began to be investigated: Earl (1961) and Meyer and Vogt (1961) 
found electrons which have proton-like energies, but which, in the BeV region, 
have only a small fraction of the proton intensity. Also, DeShong, Hildebrand, 
and Meyer (1964) later found that the electron flux is partially composed of 
positrons. In addition, Kraushaar et al. (1965) set a new upper limit to the high- 
energy gamma-ray intensity. Since one can assume that some of these electrons 
and gamma rays may be primary cosmic rays and some may be secondary to 
cosmic-ray interactions, it follows that a great deal of new information about the 
origins and behavior of cosmic rays in the galaxy might come from a study of these 
rarer components. Most recently, a component of relativistic electrons in the 
few-MeV energy region was found in interplanetary space by Cline, Ludwig, and 
McDonald (1964). If we assume, as a working hypothesis, that the electrons of 
these extremely low energies are also of cosmic origin, we may learn something 
more by comparing their properties with those of the higher-energy electrons, 
protons and other cosmic rays. It is the purpose of this paper to discuss the 
possible sources of these particles. 
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Cosmic rays are presumably created when very low-energy particles are in- 
jected into some region where they are then accelerated; if this region is not 
interstellar space itself, they may temporarily be stored in or near that region 
before propagating through the galaxy; finally they are modulated in the solar r 

environment before being detected. Certainly, several processes may compete 
for production, and others may take place in uncertain chronological order. 

The most fundamental model of interstellar acceleration, introduced by Fermi 
(1949), and the most quantitative model of solar modulation, described by Parker 
(1963), and combinations and variations of these, may possibly describe a great 
share of cosmic -ray origin and propagation. Since Fermi's model depends only 
on the particle's total energy per unit mass, which is a function of velocity 
(independent of charge or mass), and since Parker's model also depends on 
velocity alone at asymptotically low particle rigidity, we choose to represent all 
quantities in terms of the total energy per unit mass, y. 

The first measurements of the low-energy electrons are shown in Figure 1, 
in which the intensity, (dj/dy) e , is plotted versus y . These particles were 
found in interplanetary space with Explorer XVIII. The argument supporting a 
non-local origin of the electrons, put forth by the observers, rests on the fact 
that they undergo systematic intensity modulations of greater magnitude than any 

« 

possible cosmic-ray parent component. The actual intensity is uncertain, since 
some fraction of the steady or baseline intensity may be local or instrumental; 
thus, the true intensity is between the magnitude of a typical modulated increase 
and the total magnitude of baseline plus the increase. Both the steady and in- 
cremental intensities are shown in the Figure; a typical increase is apparently 
energy-independent and is about 50 percent the steady value. 
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Also shown are the electron measurements in the higher-energy region; not 
all observations known are included, but those are displayed that are the only 
data available in a given energy interval (Schmoker and Earl, 1965, Agrinier et 
al. , 1964, Daniel and Stephens, 1965), the most recent data with the best 
statistics (L'Heureux and Meyer, 1965), and the original data (Earl, 1961, Meyer 
and Vogt, 1961). It is inferred from the measurements on the positron to elec- 
tron ratio (DeShong, Hildebrand and Meyer, 1964, Hartman, Meyer and Hilde- 
brand, 1965) that approximately one third to two thirds of the electrons in the 1 
BeV region are from a source other than m -p.-e decay. Since we can expect 
that the n-fi-e differential spectrum peaks in the 20 < y < 200 region, all of the 
very low-energy electrons must be from another source. The more recently 
obtained higher-energy data can be fitted to a common curve, but it departs from 
the low-energy fit, shown in Figure 1 as a dashed line, by an amount between one 
and two orders of magnitude at the high energies, where y > 1000. Whether, 
after the v -)j.-e component has been subtracted out, the remaining spectrum is 
of one smooth form, or is the sum of two or more separate spectra, cannot be 
determined until after the differential intensities in the 10 < y < 100 and y > 104 
regions are known and after the v-n-e component itself is accurately determined 
by calculation and measurement of the e + /e~ ratio. The sources of the elec- 
trons of y £ 50 most likely pose a quite separate question from that of the 
higher-energy electrons. The very low-energy electrons are shown again in 
Figure 2a compared with a representation of the composite primary cosmic-ray 
proton intensity, observed by McDonald and Ludwig (1964), and the primary 
alpha-particle intensity, observed by Fan, Gloeckler, and Simpson (1965). 

These electron, proton and alpha measurements, made on the same satellite, 
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are shown together with balloon measurements by Balasubrahmanyan and Mc- 
Donald (1964) and by Fichtel, Guss, Kniffen, and Neelakantan (1964). The proton 
and alpha intensities fit a common curve when the alpha values are multiplied by 
a constant factor of %5, to take into account the relative abundances. This com- 
posite representation of the 1963 proton and alpha spectra, which was introduced 
recently by Balasubrahmanyan, Boldt, and Palmeira (1965) fits the data all the 
way through this low-energy region; adapted to our units, it is 

~p~ = -35 (y - l) 1-5 (y - .47)' 4 particles/cm 2 sec-sr-unit y. 

The similarity between the observed electron and proton differential velocity 
spectra has been discussed in a paper by Brunstein and Cline (1965), from which 
this figure was taken. They point out that after solar demodulation is taken into 
account by use of the model of Parker (1963) the electrons and protons fit even 
more closely the same spectrum, as shown in Figure 2b. The fact that this 
spectrum is a simple power law in total energy, which, for the protons at least, 
is continuous from a nonrelativistic y of 1. 02 up to the extreme relativistic region, 
is suggestive of a Fermi acceleration. However, the fact that this fit is accom- 
plished before the usual correction for the 2. 5 g cm' 2 galactic path length is 
made, means that either it is accidental or that electrons and protons are indeed 
propagated in a spectrally neutral fashion. As outlined by Brunstein and Cline 
(1965), several possibilities present themselves; low-energy electrons and 
protons may be Fermi-accelerated from neutral material in the galactic medium 
in such a manner as to over-compensate ionization losses, or the 2. 5 g cm' 2 
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path with its attendant spectral alteration and fragmentation probabilities may 
apply only to the heavies, or the electrons and protons may have a metagalactic 
origin, in which case during their travel to the solar system they encounter very 
little material. It is clear that these questions cannot at present be answered. 
What we will investigate in this paper is the possibility that the low-energy 
electrons are secondary to galactic cosmic rays, and the relation of this possi- 
bility to the question of a primary electron source. 

Hayakawa and Okuda (1962) originally suggested that high-energy electrons 
should arise from 7 t - /a - e decay following meson production in collisions of 
cosmic rays with the interstellar gas. Their calculations, and those of Jones 
(1963), Pollack and Fazio (1963), Ginzburg and Syrovatskii (1964), Gould and 
Burbidge (1965), and others showed that indeed a certain fraction of the observed 
electrons of y ~ 50 may originate in this manner. The high-energy positron to 
electron ratio, found by DeShong, et al. (1964) and by Hartman, et al. (1965) to 
be low, indicates that less than half of the observed electrons in the y £ 200 region 
are secondaries from this process. Brunstein (1965) showed that Coulomb col- 
lisions could similarly account for only a fraction of the observed electrons with 
5 < y < 50 unless a Fermi-like post- acceleration was incorporated. The 
majority of electrons in all observed energy regions are therefore of tentatively 
unexplained origin. In the calculations that follow we investigate in detail some 
of the various production mechanisms which may be responsible. 



II. SECONDARY GALACTIC ELECTRONS 


a. Knock-on Electrons 

The most certain process by which cosmic rays can produce low-energy 
relativistic electrons through interactions with the interstellar medium is that 
of Coulomb collisions. In this process, energy is transferred to an atomic 
electron in great excess of its binding energy, and it recoils in billiard-ball 
fashion; free plasma electrons are turned into cosmic-ray secondaries in the 
same way. These knock-on electrons are thus a necessary consequence of the 
traversal of cosmic rays through the galactic medium. 

We use the cross section for knock-on production calculated by Bhabha 
(1938). The particular form of the cross section incorporated is that for spin 
Vi , in order to obtain quantitatively accurate results to very high energies, 
y % 10 5 orlO 6 . The form used by Brunstein (1965) was for spin 0 but was 
sufficiently accurate for the 5 ~ y ~ 50 region. The differential cross section 
for the production of an electron having total energy per unit mass in dy e at y e 
by the collision of a cosmic ray of particle species j having normalized energy 
y- with a target of charge Z r and atomic number A. is 


j) d r e = 

cm 2 / gram, (1) 


r 


2ttN Z. r 2 Z. 2 

o i e j 


< 


s \y. + 


s + 1/s 


Ail 1 -/]* 2 ) t^e" 1 ) 2 ( y e — 1 ) / / 


2y 2 


dy„ 
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in which N o is Avogadro's number, r e is the classical radius of the electron: 
r e = e 2 /me 2 =2.82 x 10' 13 cm, and s = m e /m p ^ 1/1836, and for which the 
maximum transferable energy is 


y„ 


= l + 


( r i 2 - 1 )/( s { y i 


( 2 ) 


We make the usual approximation that the elemental abundances in the inter- 
stellar medium relative to hydrogen, a. , can be represented by taking that of 
helium to be about 0. 10 and by ignoring the higher-Z components. This approxi- 
mation introduces a negligible error, particularly since the cross section varies 
only as Z i /A i to the first power. The contributions by the various nuclear species 
in the cosmic radiation are more important, varying as the square of z.- Using 
Hie cosmic-ray abundances relative to protons, b . , given by Ginzburg and 
Syrovatskii (1961), we estimate that the total contribution from primaries of 
charge Z. >_2 will be about an additional 0. 75. Thus 

♦(n) = f v f b iV*ihi) ~ 1 - 75 ®iW 


The source density of electrons due to galactic knock-on production is given 
by the following integral, 

, ^ P2 d J ( r) 

Q(y e ,r) = 1. 75p (r) 4 77 (^ ’ ( 3 > 

JJ j 7 P 

electrons/cm 3 -sec-uni t 7 , 
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in which p ( r ) is the interstellar density in g/ cm 3 as a function of galactic 
position r , ^ is the differential cosmic-ray proton intensity in particles/ 

' p 

cm 2 -sec-sr-unit y , and in which the limits of integration are functions of 


electron energy y l 2 


7 


1 . 2 


K) 


We ignore any possible dependence on galactic position , r , and use the 
cosmic-ray proton spectrum discussed by Brunstein and Cline (1966), which is 
a smooth fit to the high-energy spectrum, summarized by Ginzburg and Syro- 
vatskii (1964) , to the intermediate energy spectrum' (McDonald and Webber, 
1961) and to the solar demodulated low-energy proton spectrum: 


dj/dy p - 1. 1 y p ' 2,5 protons/cm 2 -sec-sr-unit y. 


( 4 ) 


In order to evaluate the limits of integration in equation (3), we use the con- 
straint on the cross section expressed by equation (2). Since the maximum 
possible energy is unlimited, y 2 = °°, while y is determined by solving the 
inequality y e < y max for y p . The result is that y p > y x , where 

yi = ifre ~ l ) + I 1 + + s2 )( y e " !) + \{y e - O 2 } • 


It can be shown that y l is greater than the largest root of (hj (y e , y p ) as a func- 
tion of y p , which means that (hj is positive throughout the range of integration. 
Performing the integral of equation (3), we obtain 
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r 


2 77 N r 2 


Q(y e ) = 1.75(477 P )(1.1) ^ 7 _ jyi 


J 2 4 

It 7, 


1 + 


(y, - !) 


(1 + s 2 ) (y. - 1) ' 
2 


77 

~2 




7 ,’" “ 1 

y? + i 


+ arctan 


7i 





electrons/cm 3 -sec-unit 7, 
in which s and y 1 are given above. 

Since the cross section is valid over the range 1 < 7 e < 10 6 and the cosmic- 
ray spectrum (which strongly depends on the modulation model only in the region 
7 p < 2 where the electron production is of minor importance) is reliable over 
the region < 10 6 , the resulting source strength Q(/ e ) is valid over that entire 
region of . An approximation to the expression of equation (5), which allows 
for mathematical simplifications and gives an excellent numerical fit in the 1.2 
< 7 e < 10 6 region of physical interest is plotted in Figure 3. It is 

Q (y e ) = 4.9 y° (y e “ 1) 2-76 eleetrons/cm 3 -sec-unit y. (6) 
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In order to check this result we can instead integrate the cross section over 
the proton spectrum of McDonald and Webber (1961). The form they fit to the 
observations is J (>R) - 0.40 R“ 1,25 protons/cm 2 -sec-sr, in which the rigidity 
R is measured in BV. We differentiate this spectrum and change the units to 
obtain 


dj = dj dR 

d y dR dy 
p ' p 


°-54 7 p (y p 2 


1 - 1.625 


protons/cm 2 -sec-sr-unit y. 


The computed source density using this spectrum in its restricted region of 
validity 2 £ y p ~ 10 2 t is nearly the same numerically as that of equation (6). 
The source density used by Brunstein (1965) for the energy interval 
5 < y e < 50 was Q' ( E g ) = 0.91 p E e * 2,625 electrons/cm 3 -sec-MeV, which can 
be transformed: 




dE 

dy 


0.91 p (y e - 1 )* 2 ' 625 (me 2 ) 


1.625 _ 


2.7 



2.625 


electrons/cm 3 -sec-unit y. 


This source density, in its restricted region of validity, is also very close to 
that of equation (6). 
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b. Neutron-Decay Electrons 


Cosmic rays produce secondary relativistic electrons also by means of 
nuclear reactions. In this process intermediate excited nuclei or neutrons are 
produced which in turn beta-decay. This two-step process systematically shifts 
the kinetic energies involved from the BeV region to the MeV region by giving 
the decay electrons a distribution in energy per unit mass resembling that of the 
parent neutrons or excited nuclei. Since the lifetime of the neutron against de- 
cay (1000 sec. ) is infinitesimal compared to its lifetime against interaction, each 
secondary neutron will give rise to one secondary electron. 

There are at least three sorts of cosmic-ray interactions with the inter- 
stellar medium which produce secondary neutrons. (1) Pi + mesons are produced 
with neutrons in the reaction p + p- p + n + rr + . This process has the lowest 
energy threshold of any meson-producing reaction and excess tt + production 
continues to occur in the higher-energy region of multiple 7 t + and -rr~ meson 
production. Except for the case of the relatively rare process p + p - d + tt* 
which competes for t r + production, charge conservation insures that the total 
number of excess tt + mesons, integrated over all energies, is equal to the num- 
ber of beta-decay electrons at production. (2) Neutrons also are made by the 
evaporation of excited compound nuclei formed in cosmic-ray interactions. Since 
about 10% of interstellar material may be helium and about 15% of cosmic rays 
are alphas or heavier nuclei, nearly one quarter of all cosmic-ray interactions 
with the interstellar medium involve compound nuclei capable of excitation. 
Neutron emission is the most likely method of the deexcitation of compound 
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nuclei. Although the neutrons are thermal in the reference frame of the nucleus, 
the cosmic-ray motion insures a spectrum up to relativistic energies in the 
observer's frame of reference. (3) Neutrons can also be made by a nuclear 
knock-on process analogous to that for atomic knock-on electrons. Neutrons 
made in this way have somewhat greater kinetic energies than do evaporation 
neutrons since less energy goes to the remaining nucleus. 

The existence of these three processes requires that a secondary component 
of electrons should arise from nuclear reactions and also requires that this 
component should be characterized by energies in the MeV region. The 
spectrum of beta-decay electrons resulting from interstellar interactions can be 
reasonably estimated, but unlike that of knock-ons, cannot be calculated ac- 
curately since exact expressions for the nuclear reaction cross sections do not 
exist. 

The probability of a neutron at rest emitting an electron within dyj at y 
in its decay reaction n - e" + p + v is 

f i ( 7 i ) dy i - 0.614 7^7^- 53 - y 1 j 2 dy 1 elect ron/neut r on. 

This beta-decay spectrum covers the kinetic-energy region from 0 to 782 kev, 
i. e. , from a y x of 1. to 2.58. The spectrum of electrons in the observer's 
frame, f 3 (y 3 ) depends on both the spectrum of electrons relative to the parent 
neutrons f x (yj, and the spectrum of neutrons relative to the observer f ^ ( y J. 
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Analogous to the calculation by Jones (1963) for n - /x - e electron productions, 
we have, for cases in which the distribution f 3 (y 3 ) of electron energies in the 
decay frame does not depend on the spectrum f 2 [y 2 ) of decay reactions relative 
to the observer, and for cases of isotropic collisions or decays, 


f 3 W d ^3 



f i (yQ d yi c y 2 f 2 (y 2 ) d y 2 

( 2 ~ l Y ^ y 2 _ (y 2 2 ~ *)* 


(7) 


Electrons produced in dy x at y 1 are uniformly spread over the interval 


y 2 y 1 ~ [y 2 ~ 1 ) w (^i 2 ~ O’ 4 £ y £ y 2 y x + (y 2 2 _ 0* [y 2 ~ l Y A 

(Rossi, 1954). Thus, the contribution to the interval dy 3 at y 3 from dy 1 at y x 
is integrated over the portion of the spectrum f x (y x ) which contributed to dy 3 , 
using the connection between y x and y 3 expressed by the limits of integration 
on y 2 : 


y ■ 


± 


y\y^ 


± ( 7l 2 - l)'"(y 3 2 - l)’ 


A neutron spectrum f 2 (y 2 ) can be replaced by a discrete series of intensi- 
ties of monoenergetic neutrons in order to obtain a numerical result. The 
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beta-decay spectra resulting from several monoenergetic neutron distributions 

/•GC 

are illustrated in Figure 4. The area f 3 (y 3 ) dy 3 under each curve is 1. 

The decay spectrum of the electron is given by equation (5), but the neutron 
spectrum relative to the galaxy f 2 (y 2 ) depends, in turn, on both the neutron 
spectrum relative to the parent interaction f 4 (y 4 ) and the spectrum of cosmic- 
ray interaction, relative to the galaxy f g (y 5 ) . The secondary neutron spectrum 
f 4 (y 4 ) is a combination of at least three kinds of production spectra, as outlined 
earlier, and is unknown. However, it can be assumed that in the interaction, or 
center of momentum, frame, most of the secondary neutrons have energies in 
the few-MeV region, low compared with BeV-like cosmic-ray energies in the 
galactic frame. As shown by Jones (1963), a primary power-law spectrum, 
when such an inequality is obeyed, produces again a power-law spectrum of 
secondary products with the same shape but a different normalization. In the 
absence of detailed knowledge of the neutron production by cosmic rays, we can 
therefore assume that the neutron spectrum in the galactic frame reflects the 
cosmic ray spectrum with a multiplying constant containing the cross section 
for neutron production. We further simplify by using an energy-independent 
cross section of 10 millibams, based on inspection of the data tabulated by 
Pollack and Fazio (1965). This appears to approximate the cross section for 
excess 7 t + over v~ production, and, as discussed earlier, we can approximate 
by assuming that each excess tt + is accompanied by one neutron. The source 
density of secondary neutron-decay electrons due to interactions of primary 
cosmic-ray protons is thus 

^i(^3’ r )^^3 = 47ip (?) f 3 (y 3 ) dy 3 electrons/cm 3 -sec, 
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in which f 3 (y 3 ) is related to f 2 (y 2 ) by equation (7) and we approximate 


f 2 (^2) 


d J 2 

% N q cr — dy 2 electrons/sec-sr-gram. 


Here, N o - 6.02 x 10 23 per gram, a =: 1.0 x 10" 26 cm 2 electron/proton and 
dj 2 /dy 2 

is the cosmic-ray spectrum of protons/cm 2 -sec-sr-unit y given by 
equation (4). This calculated electron source density is found to be spectrally 
at least one decade below the knock-on source density. The beta-decay source 
density is thus so much less than the knock-on source density that, although its 
calculation is much less accurate, we assume its contribution can be neglected 
in the calculations that follow. 

c. Galactic Energy- Loss Mechanisms 


Electrons produced by the various cosmic-ray processes in the galaxy must 
subsequently lose energy in the galactic medium before they are observed. Since 
the medium is composed of matter, fields and radiation and has a boundary, the 
energy-loss mechanisms include ionization, bremsstrahlung, synchrotron radi- 
ation, inverse Compton effect and leakage into metagalactic space. Some of 
these are of greater importance than others; several authors have considered 
these effects for the higher-energy electrons from the 77 -» /x - e process (e. g. , 
Gould and Burbidge (1965), Ginzburg and Syrovatskii (1964), Jones (1963), 
Hayakawa and Okuda (1962) and others) and we use their commonly accepted 
values of astrophysical parameters in order to best compare results. 
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Ionization loss in the interstellar medium can be taken to be energy-inde- 


pendent over the electron energy range of interest here; although the loss does 
increase slightly at the highest energies, other loss mechanisms produce a 
much greater effect in this higher energy region. Hayakawa and Kitao (1956) 
have calculated the energy loss as a fraction of degree of ionization of the 
medium ; if we assume the galactic medium to be mostly hydrogen and to be 


about 90 per cent neutral, then we can use a value of 


dE 


d s 


5 Mev cm 2 /gram. 


The energy loss rate is directly proportional to the interstellar matter density; 


so if we then assume a constant p 


halo 


2 x 10* 26 g/cm 3 , we can write 


d y 
d t 


\ cp 

dE 

Aon m C 2 

d s 


6.0 x 10* 15 = k unit y/sec. 


(8) 


In the disk, the value is about 50 times as great. It is interesting to note that 
the rectilinear range in the halo of even a low-energy relativistic electron is ~ 
10 25 cm, much greater than the galactic diameter; in the disk where p ~ 10* 24 
g/cm 3 , it is somewhat less than the diameter. 


Bremsstrahlung is the radiation which occurs when electrons decelerate in 
matter; the loss rate is proportional to the density and, using the values of the 
constants previously assumed, we can write 

4c r e 2 Z r 2 /;N c 


d y 
dt 


brem 


137 \ 


O , / 183 ' 
lnl 

W 3 , 


(y - 1) =4.35 x 10* 18 y unit y/ sec. 
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This loss rate increases with energy, exceeding the rate due to ionization loss 
at energies above 700 MeV. 

Synchrotron radiation, which occurs when electrons accelerate moving 
through magnetic fields, is an energy-loss mechanism which is independent of 
the matter density p , but depends on the square of the galactic field component 
B ± perpendicular to the electron's velocity. If we assume the mean B ± to be 
3 x 10' 6 gauss, we can write 


d y 
d t 


2 cr e 2B ± 2 {y 2 ~ X ) 


syn 


3 me ' 


= 1.73 x 10 


- 20 


y 


by 2 unit y/ sec. (9) 


In this case the rate increases by approximately the square of the energy, ex- 
ceeding the rate due to ionization loss at energies above 300 MeV. Since the 
energy loss due to synchrotron radiation is greater than that due to bremsstrah- 
lungatall energies for which either is more important than ionization loss, we 
will ignore the effects of bremsstrahlung in our calculation. 

The energy loss due to inverse Compton collisions of electrons with inter- 
stellar photons is proportional to the radiation density and, in the region of 
interest here, to the square of the energy. If we assume that the typical galactic 
halo radiation density p r is approximately 

0.1 eV/cm 3 = 1.6 x 10~ 13 erg/cm 3 , 
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we can write 



comp 


-p. 


B 2 /4 ii 


d y 

dT 


= 3.9 


10 


- 21 


y 2 unit y/ sec. 


syn 


( 10 ) 


This loss rate is less than, but has the same energy dependence as, the synchro- 
tron rate over all the energies of interest here, 1 < y & 10 s . We shall add the 
contributions from the two effects in the work that follows. 

The escape rate due to the diffusion of the electrons into metagalactic space 
is even less certain than the loss rates considered above, since the galactic 
trapping of electrons in the energy range of interest here is unknown. Electrons 
of these very low rigidities may be trapped in the galaxy for a time similar to 
the lifetime usually assumed to be representative of the high-energy cosmic- 
ray protons, or they may be trapped for a duration orders of magnitude longer. 
High-energy electrons are usually assumed to have the same lifetime against 
escape as do cosmic rays (Hayakawa and Okuda, 1962; Ginzburg and Syrovatskii, 
1964; Gould and Burbidge, 1965), namely, T L = 3 x 10 15 sec. We will leave 
this an open question, and assume T L for these low energies to be energy inde- 
pendent and to be between 3 x 10 15 sec. and 00 . For purposes of comparison, 
an effective energy loss rate is taken to be the total energy of the particle divided 

by the lifetime: 

y 

~ 1 3.3 x 10" 16 y uni t y/ sec . 
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The loss rate increases with energy, in the limit exceeding that rate due to 
ionization loss at energies above 9 MeV, and being exceeded by that due to 
synchrotron loss above 10 BeV. It thus may dominate ionization loss and synch- 
rotron and other losses over most of the energy range of interest in this paper. 
This process, for purposes of computation, produces a particle sink, rather 
than an alteration of energy, and we therefore use 


Q l O) - 


n ( 7 ) 

T. 


< 3.3 


10 


- 16 


n = \n elect rons/cm 3 - sec-unit y 


d. Calculation of Secondary Galactic Electron Intensity 

The intensity of low-energy electrons produced in the galaxy by cosmic 
rays and subject to the energy losses outlined above is calculated by use of the 
continuity equation for the particle density in y space: 


Bn . 

+ div r ~ S 


Here n(y, t)dy is the spatial density in dy at y of electrons/cm 3 and S(y) 
d T the source production rate of electrons/cm 3 -sec. The continuity equation 
is rewritten then, as 

Bn B 4. dy i 
Bt + By n \? dt 
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in which the right hand side is the algebraic sum of the particle sources and 
sinks Qj , and 


3 

3 y 



dy d t 


is the net increase in the spatial density of electrons in d y at y due to the 

various energy loss or gain mechanisms. Inserting the production and escape 

dy. 

terms and the energy loss terms — j-j- discussed earlier, we have 


3 n 
3T 



by 2 ) 


3 n 

— + ( k ~ 2 b 'X) n 


Qx 00 


(ii) 


in which k, b and A. represent the positive numerical coefficients outlined 
above for ionization loss, synchrotron loss, and escape, respectively, and q i 
refers to the knock-on source density. The solution to this equation for n (y ) is 
shown in the Appendix to have various forms , depending on which various 
energy-loss mechanisms are taken to be dominant. Since knowledge of the 
astrophysical parameters is uncertain, we exhibit here all the solutions. 

First, in the limiting case at the lowest energies when only ionization 
losses and escape losses are relevant, 
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>Ay/k 


J -7 + tk 

Q(y') e dy electrons/cm 3 -unit y 

'V 


n e (y) 


If we use for Q(y') the knock-on source expression of equation (6), and take 
t - 00 , the integral can be analytically evaluated for A - 0, and can be numeri- 
cally evaluated for finite A. . The resulting differential intensities 



c n 


e 


477 


electrons/cm 2 -sec-sr-unit y 


are plotted in Figure 5. 

Second, in the much higher-energy region when synchrotron losses pre- 
dominate, 


d i , 


c e 


4 77 (k + by 2 ) 


f (k,A,b,y) f C 
(k + h7 2 ) J-v 


Q(y , )e -f dy‘ electrons/cm 2 -sec-sr-unit y, 


in which 

k b l/t y 

f (k, k, b, y) - — — arctan — 

(kb)* k* 
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In the event of no escape losses, we have 


d J e 

dy 


4 77 (k + by 2 ) 


These calculated intensities are also plotted in Figure 5. Comparison of the 
calculated secondary intensities with the intensity observed by Cline, et al. (1964) 
shows that the calculated absolute intensity is too low. Further, using the 
modulation model of Parker (1963), the observed value would be a factor of 2.7 
higher, increasing the discrepancy. We claim that this discrepancy is not 
trivial; i. e. , the lack of knowledge of astrophysical parameters is not a con- 
tributing factor in this lack of fit. In the limiting case in which only ionization 
losses contribute, the secondary electrons are in an equilibrium state such that 
the matter density cancels out: 


d 1 e 

dy 


b= o 
A - o 


c 

47rk 


00 

Q(y')dy‘ 

J y 


c / 4.9 p 

4 7Tk \ 1. 76 (y - l) 1 - 76 


electrons/cm 2 -sec. -sr-unit y , in which 1< is directly proportional to p. Thus, 
since the ionization losses and knock-on cross section are well known, the cal- 
culated intensity in this case depends only upon the assumption of the galactic 
universality of the cosmic-ray intensity. 
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HI. SECONDARY METAGALACTIC ELECTRONS 


a. Production Mechanisms 


Secondary electrons can be made in metagalaetic space by cosmic-ray inter- 
actions with that medium. Although the character of the intergalactic cosmic- 
ray intensity and of the medium are poorly known, we can still expect that the 
secondary low-energy electrons will be in an equilibrium similar to that in the 
galaxy. As was seen in the previous section, the fundamental quantity to know 
in order to calculate the secondaries is the primary cosmic-ray intensity. Some 
theorists assume that the metagalaetic intensity is about the same as that ob- 
served in the galaxy; for purposes of calculation we will take it to be an upper 
limit. 

Although the degree of ionization of the intergalactic medium is uncertain, 
we can make the usual assumption that it is essentially totally ionized; the cross 
section for knock-on production by Coulomb collisions with free plasma electrons 
is essentially the same as that for bound electrons. Nuclear interactions of 
cosmic rays with free protons will produce neutron-decay electrons as well, but 
again to a lesser extent; in fact, evaporation neutron production will be even less 
competitive in intergalactic space if we assume that both the cosmic rays and the 
medium contain negligible portions of alpha particles. 

We then take the metagalaetic secondary source density to be given by 
4 9 

Q m (^e) ~ i 75 P (y e ~ l) ‘ 2-76 £ 2.8 p (y e ~ l) ‘ 2-76 electrons /cm 3 - sec -unit y, 

in which p is the density in g/ cm 3 . 
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b. Metagalactic Energy- Loss Mechanisms 


Secondary electrons produced in intergalactic space lose energy because of 
their interaction with that medium, and also vanish because of the universal 
expansion. Assuming knowledge of the physical parameters involved, these 
losses can be calculated. In addition, there may be a distortion of the energy 
spectrum due to the effects of penetration into the galaxy, just as there is into 
the interplanetary region. 

Ionization losses in a fully ionized medium are replaced by losses due to 
plasma oscillations; the cosmic-ray electron loses energy by interacting with 
the plasma ions and electrons to a greater extent than it would in a neutral 
medium since the shielding of the individual charges is reduced. As shown by 
Hayakawa and Kitao (1956), for an ionized medium we have 


dE 
d s 


10 Mev -cm 2 / gram. 


|mg 


If we assume a universal proton density of 10* 5 /cm 3 , or an equivalent matter 
density p mg - 1.6 x 10* 29 g/cm 3 , using equation (8) we have for the energy 
loss rate 

“ 1.0 x 10' 17 unit y/ sec. 

ion 
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It is interesting to note that the rectilinear range of a low-energy relativistic 
electron is about 10 28 cm, close to the Hubble radius. Thus, electrons produced 
universally can reach the galaxy. The ratio of metagalactic to galactic brems- 
strahlung losses is the same as that ratio for "ionization" losses. This is because 
both interactions depend linearly on the matter density, and because a factor of 
2 is introduced in each loss when one substitutes an ionized medium for a nearly 
neutral one. Thus 



7.0 x 10* 21 y e unit y/sec. 


the bremsstrahlung loss equals that due to plasma excitation at 


= 14 


103 


or about 700 MeV. 


Synchrotron and inverse Compton losses depend on the typical metagalactic 
magnetic field strength and photon density, which are unknown. Gould and 
Burbidge (1965) adopt 10* 7 gauss, and 6. 0 x 10* 3 eV/cm 3 ~ 10* 14 erg/cm 3 
respectively, as typical assumptions. For purposes of comparison, then, we 
have from equations (9) and (10) 


dy\ 

^ — J zz 2.0 x 10* 23 y 2 unit y/sec 

' syn 
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and 



2.4 x 10‘ 22 y 2 unit y/sec. 


Inverse Compton losses thus dominate synchrotron losses at all energies. Since 
plasma-excitation losses dominate inverse Compton losses for all energies up to 
100 MeV, well above the energy region where bremsstrahlung competes, Comp- 
ton effect is the only high-energy loss mechanism we will consider. 

The magnitude of the effective particle sink represented by universal ex- 
pansion, the red-shift loss, is uncertain since the diffusion of electrons in 
intergalactic space cannot be accurately represented. We will therefore, 
analogous to the galactic case, take as an upper limit the rectilinear loss due to 
the three-dimensional expansion, characterized by a lifetime one-third the Hub- 
ble factor: T L - 1/3H. The equivalent energy loss rate is 

\ y 

J £ ~~~ £ 1.0 x 10 * 17 y unit y/sec. 

/l X L 



This loss rate is exceeded by that due to inverse Compton effect above 20 BeV 
but it actually exceeds that due to plasma excitation at all energies above y = 1. 
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F or our computation we represent this process as a sink of particles and so 
- Q L (y) ~ ^ £ 1.0 x 10' 17 • n electrons/cm 3 -sec-un.it y . 

c. Calculation of Secondary Metagalactic Electron Intensity 

The intensity of low-energy electrons produced in intergalactic space by 
cosmic rays with an assumed universal intensity is calculated by use of the 
equations derived earlier. Using the values stated above for intergalactic 
plasma-excitation and inverse Compton losses, we solve for two extreme cases, 
namely, that for the rectilinear red-shift loss and that for none. These are 
plotted in Figure 5. It is seen that the metagalactic flux, even for a primary 
cosmic-ray flux equal to that observed here, is at least an order of magnitude be- 
low the calculated galactic electron flux. The analogous result for the high-energy 
secondaries was obtained by Gould and Burbidge (1965). Thus, since this cal- 
culated metagalactic electron intensity is an upper limit, it represents a negli- 
gible addition to the galactic secondary flux, and we must look elsewhere for a 
possible source for the experimentally observed electrons. 

IV. ENERGY INPUTS TO SECONDARY ELECTRONS 

In this section we consider the possibility that the observed low-energy 
electrons are galactic secondaries which have been accelerated after their 
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creation. Two mechanisms for increasing the electron intensity to that observed 
present themselves, one in the solar system producing a local intensity increase, 
and the other producing an intensity increase throughout interstellar space. 

a. Heliocentric Electric Field 


The possibility of the presence of a substantial potential of the solar system 
with respect to the galaxy as a whole has been considered by several authors. 

In particular, attempts have been made to explain the solar modulating effects 
of primary cosmic-rays in terms of such a potential. Earlier work by Ehmert 
(1960) required a positive potential of 1 or more BV to account for the observed 
modulation. A more recent treatment by Freier and Waddington (1965) required 
50 MV to obtain quantitative agreement between their model and the observed 
spectrum. It is nevertheless difficult to envision how even such a modest 
electric field could be maintained in the presence of the interstellar gas which 
is, undoubtedly, somewhat ionized. 

In this section, we consider the effects of such a heliocentric electric field 
on the interstellar secondary electrons, which we have shown must be produced 
with an intensity not far from that observed. 

The Liouville theorem states that the differential flux of electrons 


dj 

(4 7T r 2 d r) (4 7rp 2 dp) 
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remains constant in phase space. From this, we find that 



d Ji 

is a constant. Relating the spectrum outside the solar system,- . to that 

d J 2 

observed, ^ , we have 


££2 
d y 



ymc 2 


+ qV ) 2 





1 


where qV is the product of the charge on the electrons and the accelerating 

d J 2 _ dj , 

potential. Relating the measured intensity, defined here as [y 

qV), to the calculated secondary intensity, (ymc 2 ), we find that £ 

potential of 1. 5 million volts provides an excellent fit to the data at the low 
energies. This result is shown in Figure 6 . 


me" 


The intensity increase provided by this 1 . 5 MV field is insufficient to 
simultaneously fit the observations at the high energies. The negative electrons 
in the fractional- Be V energy region, found to be in excess of those from 77-/2 
e decay by DeShong, Hildebrand, and Meyer (1964), would thus have to be due 
to yet another source. Such may or may not be the case, but a third source is 
undesirable from the point of view of minimizing hypotheses. The observations 
may yet indicate that there could be a change in the spectral form between 20 and 
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200 MeV, but it will have to be determined whether this is due to solar modulation. 
If we assume that there is only one electron source other than the v -> /j. - e pro- 
cess, we thus need to look elsewhere for the required energy input. Further, a 
50 MV field would so grossly distort the 3-MeV electron spectrum from galactic 
knock-on production as to be experimentally evident from a vast over-supply of 
low-energy electrons. 

b. Galactic Fermi Acceleration 

The acceleration mechanism of cosmic-ray protons most generally considered 
as a possible means of interstellar energy increase is the mechanism first in- 
vestigated by Fermi (1949). In this model, a particle randomly gains energy in 
the process of colliding with moving magnetic irregularities in the medium in 
which it is confined. The probability for it to have total energy between y and 
y + d y is 


d v 

dy 


(y) dy 




l/Ta 


Ta y\ + l/Ta 


dy , 


in which y 0 is the particle's value at injection, T is the mean life of the parti- 
cles, and a - v 2 /(c 2 r), in which v is the scattering center velocity and t is 
the mean time between collisions. Thus, 

j l/Ta 

d J _ Jo 'o 

dy Ta yi + i/Ta 
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in which 


Jo " 



showing that the resultant cosmic-ray spectrum should obey a power law in total 
energy. Since this result is, in fact, consistent with observations, a mechanism 
of the Fermi type gains much support, whether operative in the interstellar 
medium, or in stellar atmospheres, or in supernovae. 

If we now investigate the possibility that electrons may also be Fermi ac- 
celerated, we treat the mechanism as an energy input in the solution of the 
continuity equation for electrons in energy space. Since, in the Fermi model, 
the energy acquired by a particle of age t is y - e at in which a is defined 
above, we have 


+ 



a ’ y unit y/ sec. 


Inserting this term in the continuity equation, we have 


yjr- ( k " a y + b ^ 2 )iT7 + ( a + 2h 'y + k ) n ~ Q(y) 


( 12 ) 
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in place of equation (11). Since there now exist both positive and negative energy 
loss terms, the solutions to this equation do not directly follow from those used 
earlier. As shown in the Appendix, the forms of the solutions vary, depending 
on the choice of the parameters k, a, b and -V. In particular, the algebraic 
sign of (a 2 - 4 kb) = ^ 2 is the determining factor; i. e. , whether a is greater or 
smaller than about 2 * 10 * 17 sec' 1 , using a halo density of 2 x 10 ' 26 g cm' 3 . 

As we shall see, any value of alpha that comes close to giving a fit to the ex- 
perimental data is much greater than this; so we can, except for the postulated 
a = 0 case discussed earlier, take (a 2 - 4kb) to be positive. 

As discussed in the Appendix, care must be taken that the behavior is cor- 
rectly interpreted near the singularities which occur at energies 


7 


t 



where the rate of energy loss changes signs. The solutions for 


d i e = cr » e 

d y 4 77 


are described for the various cases of escape and synchrotron losses being or 
not being incorporated. For y < y + , which we shall see is the energy range of 
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interest here, the general solution is 



Jy 


Its limiting form for k - 0 follows directly but for b - 0 it is more convenient 
to solve the original differential equation to yield 


d i e 

dk 



r k/a 


c 

4 77k (1 - ay/k) 1 + 


Q(y) 


(i 


ay'/ k)* /a d y' 


which also has an evident limit for k - 0. These results are plotted in Figure 
7 for various values of <*. If a is large enough it is possible to match the inten- 
sity but with a poor fit to the spectral shape. 

If one attributes the primary proton beam to Fermi acceleration and the 
further assumption is made that the. parameters a and k are nearly the same 
for protons and electrons, as shown by Fermi (1949), we have for the proton 
beam 


d y 


const. y"(i +A /“) . 
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The exponent (1 + k/a) is well established to be approximately 2.5. We now must 
have k > a instead of k < a. Here /V = 1 . 5a, and the singularity in the differen- 
tial electron spectrum is removed. This solution is shown in Figure 8 for a 
variety of values of k and a, including those used formerly. In this case, toward 
the lower energies the resulting spectrum becomes independent of k and a and 
has the same slope as the observed spectrum, but, like the solution for a = 0 
which it approaches, is too low in intensity. It appears that if the electrons are 
propagated similarly to the protons, such that a Fermi acceleration with k/a - 
1.5 is present, the secondaries from interstellar knock-on production cannot 
account for the observations. Of course, it is not obvious that either k or a 
would have the same values for electrons as they have for protons, or that K/a 
must necessarily be equal to 1.5 The electrons have a much lower rigidity than 
the protons do at a given y; they would therefore have a greater lifetime against 
escape and possibly a much different acceleration. 

We can conclude that the intensity and shape of the observed low-energy 
electron spectrum cannot simultaneously be accounted for by galactic knock-on 
secondaries. 

V. CONCLUSION 

The secondary electron spectrum to be expected from n - ^ - e decay as 
calculated by Hayakawa and Okuda (1962), Jones (1963), and others, coupled with 
the positron-electron measurements by DeShong, et al. (1964) and by Hartman, 
et al. (1965), indicate that meson-decay electrons cannot adequately explain the 
observed properties of high-energy cosmic-ray electrons. 


34 



Two of the authors have in a previous paper (Brunstein and Cline [1966] ) 
discussed the possibility that the electrons may be of a primary nature, being 
accelerated in the same process along with cosmic-ray protons. We have here 
examined the possibility that cosmic-ray electrons are secondary to the nuclear 
beam, arising from both knock-on collisions and from the beta-decay of the 
resultant neutrons from nuclear interactions in the interstellar gas and in 
metagalactic space. We conclude that neutron-decay electrons are well below 
the intensity to be expected from knock-on collisions and thus may be neglected 
as a potential source. However, at least at low energies, an intrinsically accurate 
calculation of the knock-on electron intensity shows it to be well above that 
expected from v - /x - e decay but substantially below the measured intensity 
of low-energy electrons. 

We have therefore postulated two prospective models to account for the 
observed electrons as being due to these knock-on electrons: (1) the knock-on 
electrons are further accelerated after their injection by a heliocentric electric 
field, (2) the knock-on electrons are accelerated after injection by a general 
galactic Fermi mechanism. 

The first possibility, the heliocentric electric field, requires only a modest 
1.5 MV to boost the calculated knock ons to the observations; however, it does 
not appreciably affect the higher-energy electrons where we still must face the 
question of the excessive negative electrons. Moreover, it is difficult to envision 
how even a modest potential of 1.5 MV could be maintained by the solar system 
(Lam and Sandri [ 1965] ). 

The second possibility, that of Fermi acceleration after secondary produc- 
tion, is seen to have limitations. The value of the parameter a, defined by 
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^ = ary, is varied in the neighborhood of that value which produces an 
approximate fit to the data, a £ 10“ 15 sec' 1 . When a lifetime against escape 
of about 10 8 years is assumed ( ^ = 3.3 10“ 16 sec' 1 ) the observed intensity can 
be roughly fitted, although the fit to the slope in the low-energy region is poor. 

In the high-energy region, if we assume that about half of the electrons are 
accounted for by rr-fx-e decay, we find that the remainder can also be roughly 
fitted. However, in this case, the value of a is quite large, > 10“ 15 sec., and 
is about the same as that value shown by Brunstein and Cline (1966) to describe 
the electron intensity in terms of Fermi-accelerated primary thermal electrons. 
If these parameters accurately represent galactic conditions, then we may draw 
the conclusion that the cosmic-ray beam will be composed of both secondary 
knock-on electrons and primary electrons Fermi-accelerated in essentially 
equivalent fractions. If the actual value of the parameter a is smaller, the 
secondary electrons account for a smaller fraction of the observed intensity. 

On the other hand, if the value of the ratio A./a is a constant 1.5 as is usually 
assumed for protons, then in the low-energy region a fit to the observations, 
independent of a, is impossible. 

We can conclude, in any case, that the secondary hypothesis alone cannot 
account for the observations and that at least half of the observed low-energy 
electrons have a primary origin. We note that the observations of the charge 
ratio led to a similar conclusion for the high-energy electrons (DeShong et al., 
1964, Hartman, et al., 1965). We feel, therefore, that a primary cosmic-ray 
beam of electrons is present throughout the energy region fromy £ 1 to y >10 5 , 
having about the same velocity spectrum as the cosmic-ray protons in the 
1 < y < 10 2 region. However, further observations in the y < 5 and y > 10 s 
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region are desirable to better understand the properties of these primaries, 
and more details in the 20 < 7 < 200 region are desirable to map the transition 
between the knock-on and the secondary contributions. 
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APPENDIX 


The equation to be solved is 


nH k - 


ay 


by 2 ) 


>n 


'dy 


+ (a - 2 by + A.) n 


Q(r) . 


(i) 


in which n - n(y, t) . The initial condition assumed is n(y, t - 0) - 0,and 
the solution is obtained by the method of characteristics. 


Before proceeding to the actual solution we consider two possibilities : 


1. Particles produced with certain energies lose some of their energy 

when propagating through the medium, while others, produced with energies in 

dy 

a different range, gain energy from the medium. This implies that ^7 changes 
sign, and for this to be so it is necessary that 


v 2 = a 2 - 4kb 


( 2 ) 


be greater than zero. The changes in sign occur at the energies 
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a. ± v 

2b ’ 


( 3 ) 
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When 4 kb/a 2 « 1, we have the approximations 



< 

y 




a 

"b 

k 


> 


a 


J 


2. Particles produced with any energy never gain energy from the medium. 

d y 

This implies that for all energies < 0, and hence that 


M 2 = 4 k b ~ a 2 


(4) 


be greater than zero. 


We treat these two cases separately. 


a. Predominance of Energy Increase 


Here a 2 - 4kb > 0. 


The solution for this case is 


n (y > t) 


(y-r T / -~ 1 

(y ~ y-) A/ " * 1 


/.e <y, t) 

Q(r') 


•'r 


( y - y~ 
\y - y + 



dy' , 


(5) 
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where 


g(y, t) 


y+ (7 - y~) e vt -\y~ (7 - 7 ) 
(7 ~ y~) e~ vt 


( 6 ) 


At this stage we want to investigate the behavior of N (7, t) at the exceptional 
points 7 ~ y i . One finds, by using the L’Hopital rule and the form (6) of 
f (7, t ) , that 


lim n (7, t) 

7 ^ 7 “ 


Q(7~) 


e ~(A + v)t:) < 


( 7 ) 


Therefore at y = y~ , there is no singularity in n (7, t) for all times. 

To find the limiting value of n (7, t) when y - 7 + , it is necessary to go to 
the original equation (1) and proceed there to the limit, since precisely at this 
energy there is a singularity also in the integrand of ( 5 ) in addition to the one 
outside the integral. Finally we obtain 

lim n (y, t ) = Q,(zl) /, . -(A--),' (8 ' 

r - 7* x - * ' ' ’ 
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Equation (8) shows that two possibilities may arise: 


i) k - v > 0: 


Here there is no singularity at y - y + for all times, and 


n(y + ,t 


00 ) 


Q (r + ) 

k - V 


( 9 ) 


ii) k - v < 0: 


Here there is no singularity for finite times but when t - 00 » then 


n(y + ,t 



00 # 


( 10 ) 


This situation persists when k - 0. 


In order to determine, (for case ii),the behavior of n (y, t) when a steady 
state prevails, we must investigate the rate of energy loss, 

= - (k ~ ay + by 2 ) (If) 
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and find the time intervals needed by particles produced at certain energies to 
achieve different energies. One finds that for all y < y + , a steady state will be 
reached when t - 00 ; while for y > y + , the steady state is attained after a 
finite time interval t max given by 


t 

max 


r c° 

d y 

k - ay' + b(y') 2 


1 

— In 

v 


y ~ 7 

y - y + 


( 12 ) 


(18) 


One finds that when t ^ 00 , f (y, t) - y ; and when t - t max , f (y, t ) - 00 . 
Therefore the steady state solution n s (y) is: 


n.OO 


k - ay + b y l \ y - y 


+ \A/v 

7 ” 7 
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Q(y') l 2 ^ 

y ~ y 


A /i 
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(19) 


k - ay + by 2 \ 7 “ 7 


7-7 


+ \A/v 


Q(y> (— y - 

\ y - y 


- \\/v 


dy', for y < y + 
(13) 

dy’, for y > y + 




Note that for k > v, n (y) has no singularities; while for h ± v, the only 
singularity is located at y = y* , corresponding to a pile-up of particles on the 
energy axis where synchrotron loss equals Fermi-acceleration gain. 
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The integral appearing in (13) cannot be evaluated in closed form for arbi- 


trary values of A./V . For the special case A. - 0 we have 
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lim n s (y) 
\ - 0 


n s<» 
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k - ay + by 2 


•y' ~ k/a 

Q (y 1 ) dy', for y < y + 




(14) 


k “ ay + by 2 


Q(y')dy\ for y > y + . 




If b - 0, the equation for n (y , t ) is 


Bn Bn 

- (k - ay) + (a + k) n - Q(y) 
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and the solution is 
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where 


<7,0 ■ 


There is no singularity; in fact at the point y = — = y 


o > 


we obtain 


lim n (y, t) 

7 - 7 0 
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a + A. 


Q(7 0 ) (1 - e 


at 


)• 


The steady state solution is obtained when we let t - 00 : 


n s O') 


k(l 
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ay/k ) 1 + A/a 
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for y < y 0 , and 
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n (r) 


k ( a y/k - l) 1+A/a 


Q(y') 


( , \ A/a 

c 


dy' 


r 0 = k/a 


for y > 7 0 . 


For \ - 0 (or k « a ), the steady state solution is: 


n (7) 



* k / cl 

Q(7')d y 




( 20 ) 


which is the special case evaluated by Brunstein (1965). 


b. Predominance of Energy Decrease 


Here n . 2 = 4kb - a 2 > 0. 


The solution for this case is 


n (7, t) 
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( 21 ) 
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where 


a _ / jJ- \ _ a ~ 2by ~ /x tan (/xt/2) 
2b \2b/ /j. + (a ~ 2by) tan (/xt/2) 


( 22 ) 


Equation (21) implies an electron density periodic in time, a result which 
obviously is inadmissible if not interpreted correctly. We show that beyond a 
certain characteristic time t max the expression (21) for n (y, t) loses all 
meaning. This follows immediately from equation (11), when we note that for 
the case under consideration 


dy 

dt 




< 0 


(23) 


The particles therefore always decelerate. We define t max as the time needed 
for a particle of infinite initial energy to decelerate to energy y: 



( 24 ) 
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or 


t 

max 


(7) 


77 2 2 by ~ a 

— arctan 

MM M 


(25) 


It can be shown that at t = t max , n(y, t) reaches its first maximum. 
This follows from the fact that Q(y) is given by a negative power law and from 
the result that 


1 im h (y, t) 

t - t 


max 


00 . 


(26) 


This result would not be true for a production rate with a positive power law 
form. 

Beyond t n(y, t) decreases toward zero. We therefore adopt the 
solution attained at t - t max as the appropriate steady state solution. Thus 


n s(T) 


2 A. 2 by - a ’ 

exp \ arctan 






k “ ay + by : 


(~2k 2 b y'-a 

Q(y') exp I arctan 






dy' , 

(27) 
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Again the integral in equation (27) cannot be evaluated in a closed form for arbi- 
trary K. For the case A. - 0 (or k « // ), we can write 


lim n s (y) 

\ - 0 


= n.(y) 


k - ay + by 2 


f 


Q(y' )dy'. 


(28) 


Finally, the solution for a - 0 follows from equation (27), or from (28) if 
A- - 0, unless b = 0. In that case we again solve equation (1) and find 


n(7) 


exp (Ay/k) 


b - 

a - 


0 

0 


k 



(29) 
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FIGURE CAPTIONS 


Figure 1: Observed cosmic-ray electron intensities in the 1 < y < 10 6 

region. Most of the differential intensities shown are taken directly from the 
references (Cline, Ludwig and McDonald, 1964, in which the diamonds represent 
steady values and the crosses time variations; Schmoker and Earl, 1965; Meyer 
and Vogt, 1961; Earl, 1961; L'Heureux and Meyer, 1965). The integral intensi- 
ties of Agrinier, et al. (1964) and Daniel and Stephens (1965) were converted to 
differential intensities by assuming a slope of -1.5 and an infinitely high cutoff. 

Figure 2a: The low-energy electron observations plotted, for comparison, 

with a smooth fit to a composite of the low-energy proton and alpha data, in which 
the alpha intensities have been multiplied by 5. The dashed curve is an extrapo- 
lation, beyond the measurements, of the proton and alpha fit, devised by 
Balasubrahmanyan, Boldt and Palmeira (1965). 

Figure 2b : The same data, for which the proton curve has been demodu- 

lated by the factor exp (K/ 6 ) with three values of K and for which the electron 
data have been altered with K = 1. Higher-energy proton fits (McDonald and 
Webber, 1961; Ginsburg and Syrovatskii, 1964) are shown for comparison. These 
curves were taken from Brunstein and Cline (1965). 

Figure 3 : The source strength for galactic knock-on production, Equa- 

tion 5, plotted for comparison with the approximation of Equation 6. The approxi- 
mation is used in all the calculations. 

Figure 4: Beta-decay spectra of monoenergetic neutrons. Cosmic-ray 

neutrons of characteristically BeV energies are seen to produce electrons having 
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a spectrum in the MeV region; the source strength from galactic secondary 
neutrons is found, however, to have a significantly lower intensity than does that 
from knock-on production. 

Figure 5a ; Calculated differential intensities of electrons from galactic 

secondaries are compared with sample observations. The results are shown 
with and without losses due to synchrotron radiation and losses due to escape 
from the galaxy, using a halo value of the matter density. Also shown is the 
result in which the disk value of the density was used, and losses were incor- 
porated. When no losses are used, the result is independent of the density. 

Figure 5b : Upper limits to the calculated differential intensities of elec- 

trons from metagalactic secondary production. Again, synchrotron losses and 
red-shift escape losses are varied for comparison. In all cases, the metagalactic 
intensities are below the galactic intensities. 

Figure 6 : Calculated differential intensity of electrons resulting from 

galactic knock-on production and accelerated by a heliocentric field of 1.5 MV. 

If such a field exists, the low-energy observations can be accounted for, but the 
high-energy electrons remain in excess of secondary meson production. 

Figure 7 : Calculated differential intensities of secondary electrons which 

have been accelerated by a galactic Fermi process. Various values of a are 
incorporated for which the results are plotted in a) without, and in b) with, escape 
losses. The galactic halo value of the matter density is used throughout. It is 
seen that the observations, samples of which are shown, can be roughly fitted 
with a sufficiently large value of a . . For certain combinations of the galactic 
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parameters p, b, k and a, peaks in the intensity would occur at normalized 
energies equal to k/a . 

Figure 8: Calculated differential intensities of secondaries for a fixed 

ratio of k/a = 1.5. Here the singularity disappears at the high energies, but 
the low-energy fit becomes independent of k and a and is seen to be inadequate 
to account for the observations. 
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